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ABSTRACT 



N-body simulations have shown that the dynamical decay of the young (~ 1 Myr) Orion Nebula cluster could be responsible for the 
loss of at least half of its initial content of OB stars. This result suggests that other young stellar systems could also lose a significant 
fraction of their massive stars at the very beginning of their evolution. To confirm this expectation, we used the Mid-Infrared Galactic 
Plane Survey (completed by the Midcourse Space Experiment satellite) to search for bow shocks around a number of young (S 
several Myr) clusters and OB associations. We discovered dozens of bow shocks generated by OB stars running away from these 
stellar systems, supporting the idea of significant dynamical loss of OB stars. In this paper, we report the discovery of three bow 
shocks produced by O-type stars ejected from the open cluster NGC6611 (Ml 6). One of the bow shocks is associated with the 
09.5Iab star HD165319, which was suggested to be one of "the best examples for isolated Galactic high-mass star formation" (de 
Wit et al. 2005). Possible implications of our results for the origin of field OB stars are discussed. 

Key words. Stars: individual: BD-14° 5040 - stars: individual: HD 165319 - open clusters and associations: general - open clusters 
and associations: individual: NGC 661 1 - ISM: individual objects: RCW 158 
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1. Introduction 

Dynamical three- and four-body encounters between stars in the 
dense cores of young star clusters impart to some stars velocities 
exceeding the escape velocity from the cluster's potential well 
(Poveda et al. ,1967; Leonard & Duncan 1990). The concentra- 
tion of massive stars towards the center of young clusters (e.g. 
Hillenbrand & Hartmann [19981 Gouliermis et al. l2004l Chen et 
al. p007 ) implies that the dynamical evolution of cluster cores 
is dominated by massive stars and that the majority of escaping 
stars should be massive. The escaping massive stars constitute 
a population of OB field stars (Gies 119871 ). The escape veloc- 
ity ranges from several km s"' (for loose and low-mass clusters) 
to several tens of kms"' (for compact and massive ones), so 
that the typical space velocity of OB field stars is ~ lOkms"' 
(e.g. Gies 1987). The stars with peculiar (transverse) velocities 
exceeding 30kms ' are called runaway stars (Blaauw |1961| ). 
It should be realized, however, that any star unbound from the 
parent cluster should be considered as a 'runaway', indepen- 
dently of its peculiar velocity (cf. de Wit et al. 2005 ). Although 
the peculiar velocity of some runaway OB stars could be as 
high as several 100 km s ' (e.g. Ramspeck et al. l2001l l or even 
~ 500 - 700 km s ' (Edelmann et al. l2005l Heber et al. l2008b . 
the majority have velocities of several tens of km s ' . Therefore, 
most runaway OB field stars should be located not far from their 
birth clusters. 

A (massive) star moving supersonically with respect to the 
ambient medium generates a bow shock (Baranov et al. 119711 
Van Buren & McCrav il988) visible in infrared (Van Buren et 
al. [19951 Norieg a-Cres po et al. [T997l Fran ce eta l. |2007] l and/or 
H„ (Kaper et al. 119971 Brown & Bomans [2005]) . Observations 



of runaway OB stars showed that only a small fraction (g 20 
per cent) produce bow shocks (Van Buren et al. 119951 Noriega- 
Crespo et al. 1997 ; Huthoff & Kaper 2002 ). The main reason for 
this is that many runaway stars move through the hot interstellar 
gas and their peculiar velocities are lower than the sound speed 
in the ambient medium (Huthoff & Kaper '2002). The geometry 
of a bow shock allows us to infer the direction of motion of the 
associated star and thereby trace its trajectory backwards to the 
parent cluster in those cases in which the proper motion of the 
star has not been measured directly. 



e-mail: bomansOastro . rub . de 
e-mail: vgvaram@mx . iki .rssi.ru 



N-body simulations by Pflamm-Altenburg & Kroupa (120061) 
showed that the dynamical decay of the young (~ 1 Myr) Orion 
Nebula cluster could be responsible for the loss of at least half of 
its initial content of OB stars (see also Kroupa l2004l cf. Aarseth 
& Hills 1 1972) . This result suggests that other young stellar sys- 
tems could also lose a significant fraction of their massive stars at 
the very beginning of their evolution. We therefore expect to ob- 
serve numerous signatures of interaction between ejected stars 
and the dense environment of young star clusters. To confirm 
this expectation, we searched for bow shocks around a num- 
ber of young several Myr) clusters and OB associations us- 
ing the Mid-Infrared Galactic Plane Survey. This survey, car- 
ried out with the Spatial Infrared Imaging Telescope onboard the 
Midcourse Space Experiment (MSX) satellite (Price et al. |2001| l, 
covers the entire Galactic plane within \b\ < +5° and provides 
images at 18" resolution in four mid-infrared spectral bands cen- 
tered at 8.3 yum (band A), 12.1 pm (band C), 14.7 pm (band D), 
and 21.3 pm (band E). To search for possible optical counter- 
parts to the detected bow shocks, we used the SuperCOSMOS 
H-alpha Survey (SHS; Parker et al. [2005l l and the Digitized Sky 
Survey (McLean et al. 2000). Our choice of young stellar sys- 
tems implies that the majority of ejected stars should not have 
traveled far from their birth places (so that the parent clusters 
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could be identified with more confidence) and that supernova 
explosions in massive binaries do not significantly contribute to 
the production of runaway stars (Blaauw 1961). 

We discovered dozens of bow shocks. Three are discussed 
in this paper, while the results of study of other objects are pre- 
sented elsewhere. The geometry of the three bow shocks sug- 
gests that they are driven by stars expelled from the young open 
star cluster NGC 6611. In Sect. 2, we review the relevant data 
for this cluster. The results of the search for bow shocks around 
NGC 66 11 are presented in Sect. 3 and discussed in Sect. 4. 
Conclusions are summarized in Sect. 5. 

2. NGC 6611 

NGC 66 11 is part of the Ser OB 1 association in the Sagittarius 
spiral arm. It is embedded in the molecular cloud W37 and is 
responsible for ionizing the Eagle Nebula (Ml 6), famous for its 
"elephant trunks" and "EGGs" (evaporating gaseous globules; 
Hester et al. |19961 l. Numerous distance estimates for NGC 661 1 
converge to a figure d ~ 1.75 - 2.00 kpc (Hillenbrand et al. 
[1993 ; Loktin & Beshenov |2003l Guarcello et al.|2Q07l Wolff et 
al. |2007). The age of the cluster is less well constrained. The 
presence of four early O-type stars in the cluster (e.g. Dufton et 
al. |2006) implies that its age is ^ 2 - 3 Myr, while the asso- 
ciation with the cluster of a B2.51 star (BD-13°4912) suggests 
that star formation in this part of Ser OB 1 started ~ 6 Myr ago 
(H illenbr and et al. [19931 de Winter et al. [19971 Kharchenko et 
al. 120051) . 

There is a growing body of evidence, however, that clustered 
star formation is a rapid process, so the age spread in the cluster 
should not exceed ~ 1 Myr (e.g. Elmegreen 120001 Ballesteros- 
Paredes & Hartmann 120071 ) . The evolutionary status of the ma- 
jority of OB stars in NGC 661 1 is consistent with an age of ~ 4 
Myr and in the following we adopt this figure to be the age of 
the cluster. The observed age spread could be understood if the 
most massive (and therefore the youngest) stars in the cluster 
are the blue stragglers (the rejuvenated stars formed via merg- 
ing of the less massive stars in the course of binary-single and 
binary-binary encounters; e.g. Leonard [l995l Portegies Zwart et 
al. Q999; Gvaramadze & Bomans 120081 ), while the older stars 
are the former field stars trapped in the potential well of the con- 
tracting pre-cluster cloud (Pflamm-Altenburg & Kroupa l2007l l. 
Numerous young (~ 0.1 - 3 Myr) pre-main-sequence stars sur- 
rounding NGC 661 1 (e.g. Guarcello et al. 2007 ) could represent 
the second generation of stars, whose formation was triggered 
inside the remainder of the natal molecular cloud by the activity 
of massive stars in the cluster. 

Radial-velocity measurements for OB stars in NGC 661 1 in- 
dicated that ^ 30-50 per cent are spectroscopic binaries (Bosch 
et al. I1999I I. The high percentage of massive binaries is one of 
the necessary conditions for production (e.g. via binary-binary 
encounters; Mikkola 119831 Leonard & Duncan 119901 ) of run- 
away OB stars and merged (rejuvenated) massive stars (either 
bound or unbound to the cluster; Leonard 1995 ). Two other nec- 
essary conditions are the presence in the cluster of a sufficiently 
high number of massive binaries (i.e. the cluster should be mas- 
sive) and the initially compact structure of the cluster It is likely 
that both conditions were fulfilled in NGC 6611, at least at the 
beginning of its dynamical evolution. 

A mass estimate for NGC 66 11 suggests that the cluster 
could be as massive as ^ 2x10"* Mq (Wolff et al.'2007^. This esti- 
mate was derived by assuming that ~ 160 cluster members with 
masses ranging from 6 to 12 Mq constitute ~ 5 per cent of the 
total cluster mass (for a Salpeter initial mass function). Adopting 



this mass estimate, we infer that the cluster should contain ~ 100 
OB stars, which exceeds the number of observed stars by a factor 
of 2.5. It is possible that the Salpeter initial mass function is not 
applicable for NGC 6611. Another possibility is that the cluster 
has already lost most of its massive stars by dynamical ejection 
(see also Sect. 4). 

It is likely that NGC 66 1 1 was initially more centrally con- 
centrated and that its current size [the radius of the cluster is 
- 1 - 8 pc for = 2 kpc (Belikov et al. [20001 Bonato et al. 
2006^] is the result of dynamical decay of the cluster's core and 
overall cluster expansion caused by gas expulsion due to super- 
nova explosions and stellar winds (e.g. Boily & Kroupa 2003J). 
The initially compact configuration would be consistent with the 
observation that the characteristic radius of young clusters is ^ 1 
pc, which is independent of their mass (Kroupa & Boily '2002]). 
For this radius and the above mass estimate, the escape velocity 
from the cluster is inferred to have been ~ lOkms"'. 

Proper-motion measurements for NGC 66 11 (see Table 1) 
suggest that the cluster is moving in the west-east direction 
with a peculiar (tangential) velocity of ^ 10 - 15kms"'. To 
convert the observed proper motion into the true tangential ve- 
locity of the cluster, we used the Galactic constants Rq = S 
kpc and O p = 20 km s ' (e.g. Reid 1993^ Kalirai et al. [20041 
Avedisova 120051 ). the solar peculiar motion (Uq, Vq, Wq) = 
(10.00, 5.25, 7.17) kms ' (Dehnen & Binnev [T998l ), and adopted 
a distance of J = 2 kpc. 

3. Results of the search 

3.1. Bow shocks 

We searched for bow shocks around NGC 6611 in a 12° wide 
area elongated along the Galactic plane and centered on the 
cluster's longitude (/ - 16°95). Along the Galactic latitude, the 
search was limited by the MSX coverage (see Sect. 1). At a dis- 
tance of 1.75 - 2.00 kpc, 1° corresponds to ^ 30 - 34 pc, so 
that potentially we were able to detect bow shocks produced 
by stars leaving the cluster at the very beginning of its evolu- 
tion (i.e. ~ 4 Myr ago) with a peculiar (tangential) velocity of 
$ 25 - 50kms"' . We realize that some runaway OB stars were 
ejected at higher velocities and therefore could be beyond the 
area covered by our search. On the other hand, the lower the 
peculiar velocity, the lower the percentage of stars moving su- 
personically through the ambient medium. Thus, we expect that 
only a small fraction of runaway stars located within several de- 
grees from the cluster are able to generate the bow shocks (see 
also Sect. 4) 

The visual inspection of MSX images revealed several bow- 
shock candidates. Three of them (indicated in Fig. 1 by circles) 
have a clear arc-like structure that opens towards NGC 66 11 
(Fig. 2-4), which suggests that these structures are generated by 
stars expelled from the cluster The bow shocks 1 and 2 are most 
prominent at 21.3 fim (Fig. 2 and Fig. 3), while the bow shock 3 
is visible only at 8.3 fim and 21.3 fim (more clearly at 8.3 ;um; 
Fig. 4). The bow shocks 2 and 3 have obvious optical counter- 
parts in the SHS (see Fig. 3 and Fig. 4), which are the parts of 
Hii regions, respectively, G 016.9-0 1.1 (Kuchar & Clark 1997) 
and RCW 158 (Rodgers et al. 1 19601 see also Sect. 4). 

3.2. Associated stars 

The SIMBAD database provides the spectral types of two of the 
three stars associated with the bow shocks. We found that the 
bow shock 2 is driven by the star BD-14° 5040 [with a photo- 
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Fig.l. A 6° X 6° 21.3 /urn (MSXhand E) image of NGC 6611 and its environments [including two open clusters NGC 6604 and 
NGC 6618 (M17)], with the positions of three bow shocks marked by circles. The arrows labeled as L&B and Kh show the direction 
of cluster's peculiar motion, based on the proper motion measurements, respectively, by Loktin & Bechenov ( 120031 ) and Kharchenko 
et al. ( I2005I I. The origin of the arrows corresponds to possible birth places of the cluster ~ 4 Myr ago. The thin lines represent the 
trajectories of ejected stars, as suggested by the geometry of their bow shocks and proper motion measurements (see text for details). 
North is up and east to the left. 

Table 1. Proper-motion measurements for NGC 66 11 and the three ejected stars. Two measurements are given for each object to 
indicate the uncertainties in the measurements. For each data set, the pecuhar (transverse) velocities (in Galactic coordinates) were 
calculated and added to the table. 



Object 


fi„ cos 6 


J^6 


Sources for 


Vl 


Vb 




mas yr"' 


mas yr"' 


proper motions 


km s"' 


km s"' 


NGC 6611 


0.6 + 0.1 


-0.3 + 0.1 


Loktin & Beshenov 1200311 


9.7 + 0.9 


0.7 + 0.9 


NGC 66 11 


1.6 + 0.3 


-0.4 + 0.5 


Kharchenko et al. ( 2005') 


12.9 + 4.3 


-7.7 + 3.5 


star 1 


0±7 


12 + 3 


USNO-Bl.O (Monet et al. 2003 ) 


108.7 + 39.8 


61.5 + 59.7 


star 1 


-4.3 + 8.0 


-0.9 + 8.0 


NOMAD (Zacharias et al. 2004a ) 


-18.2 + 36.9 


39.2 + 36.9 


BD-14° 5040 


5.5 ± 1.3 


-3.0+14 


NOMAD (Zacharias et al.,2004ai) 


7.5 + 13.1 


-51.7 + 12.5 


BD-14° 5040 


7.7 ± 1.6 


-4.6+1.8 


Kharchenko 1200 lb 


4.0+16.8 


-76.5 + 15.9 


HD 165319 


-1.2+1.0 


-1.4 + 0.9 


UCAC2 (Zacharias et al.|2004b') 


-9.7 + 8.7 


11.0 + 9.3 


HD 165319 


0.3 + 1.0 


-1.8 + 0.6 


NOMAD (Zacharias et al._2004a) 


-6.1+6.7 


-3.3 + 8.7 



Fig. 2. - Left: MSX 21.3 fim image of the bow shock 1 . The po- 
sition of the associated star (star 1) is marked by a circle. Right: 
SHS image of the same field with the MSX 21.3 fim image over- 
layed in black contours. Field size is 8' by 8f5, the images are 
oriented with Galactic 11 increasing to the left and Galactic b in- 
creasing upward. 



Fig. 3. - Left: MSX 21.3 fim image of the bow shock 2. The po- 
sition of the associated star BD-14° 5040 is marked by a circle. 
Right: SHS image of the same field. Field size and orientation of 
this figure are the same as for Fig|2l 



metrically estimated spectral type of 08V; Kilkenny ( 119931 )1 and 
the bow shock 3 is produced by the star HD 1653 19 [whose spec- 
tral type of 09.5Iab was derived spectroscopically by Crampton 
& Fisher dl 9741 )1. The bow shock 1 is apparently generated 
by a star with coordinates: a2000 = 18*'15™23'*.97, ^2000 - 
-13°19'35'.'8 (hereafter, star 1; see Fig. 2). This star has the vi- 
sual magnitude of 1 1.81 (McLean et al. 2000) and 2MASS mag- 
nitudes (J,H,KJ=(8.096,7.684,7.396) (Skrutskie et aL ,2006) . 
Using the extinction law from Rieke & Lebofsky ( 119851 ) and 
the photometric calibration of optical and infrared magnitudes 
for Galactic O stars by Martins & Plez (2006 ), we derived ex- 
tinction towards the star 1 of Ay - 5.6 mag and estimated the 
spectral type of this star to be 09.5III/05V (for li = 1.75 kpc) 
or 07.5III/04V (d = 2.00 kpc). The spectral type of 07.5III is 
more consistent with the adopted age of the cluster of 4 Myr 
Using the same calibration, we re-estimated the spectral type of 
BD-14° 5040 to be 07V (d = 1.75 kpc) or 06V (d = 2.00 kpc) 
and derived extinction of Ay - 3.8 mag. For both distances, the 
isochronal age of the star is less than 4 Myr, so that this star 
could be a blue straggler ejected from the cluster (e.g. Leonard 
119951 Gvaramadze & Bomans 2008). Also, we estimated the 
photometric distance to HD 165319 to be li ^ 2.00 kpc, in good 
agreement with the distance to NGC 661 1. 



Fig. 4. - Left: MSX 8.3 fim image of the bow shock 3 and the 
associated star HD 165319. Right: SHS image of the same field. 
Field size and orientation of this figure are the same as for Fig|2l 



3.3. Proper motions 

Figures 1, 2 and 4 show that the symmetry axes of the bow 
shocks 1 and 3 are misaligned with respect to the lines drawn 
through the associated stars and the center of the cluster This 
misalignment could be understood if one takes into account 
the eastward motion of the cluster (see Sect. 2 and Fig. 1; cf. 
Hoogerwerf et al. |2001| ). Assuming that there is no bulk stream- 
ing motions of the ambient gas, we found that the trajectories of 
the star 1 and HD 1653 19 (suggested by the morphology of their 
bow shocks) intersect with the trajectory of the cluster ^ 1.6 
Myr ago, if one uses the cluster proper motion measurements 
by Loktin & Bechenov (120031 ). or ~ 0.9 Myr (star 1) and 1.8 
Myr (HD 165319) ago for the proper motion measurements by 
Kharchenko et al. (2005 ). In both cases, the stars were ejected 
more than 2 Myr after the cluster formation. This inference how- 
ever would be incorrect if the ambient medium had a significant 
peculiar velocity, such that the symmetry axes of the bow shocks 
did not coincide with the direction of peculiar motion of the as- 
sociated stars. 

In Table 1 , we provide two representative proper motion 
measurements for each star, found with the VizieR Catalogue 
Service. One can see that only BD-14°5040 has a significant 
proper motion (i.e. the measurement uncertainties are far less 
than the measurements themselves). Both measurements for this 
star suggest that it was ejected ~ 1 Myr ago with a velocity typ- 
ical of classical runaway stars. It is curious that the symmetry 
axis of the bow shock produced by this star points almost ex- 
actly towards the cluster If however the star was ejected from 
NGC 66 11 and the proper motion measurements are correct, 
then the trajectory of the star should intersect with the cluster 
^ 0.30° - 0.45° to the west of its present position (for conve- 
nience, we show in Fig. 1 only one trajectory of the star). The 
"incorrect" orientation of the bow shock, in this case, could be 
attributed to the effect of peculiar motion of the ambient gas. 
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The proper motion measurements for the star 1 and 
HD 165319 are insignificant. Taken at face value, they sug- 
gest that the star 1 is a runaway in the classical sense and that 
it was able to reach its present position even if ejected from 
NGC 66 11 only ~ 0.9 Myr ago (see above). The peculiar ve- 
locity of HD 165319 (see Table 1), however, is far smaller than 
the velocity (of ~ 45 - 60kms"') inferred from the separation 
of this star from the possible trajectories of the cluster and the 
ejection time implied by the morphology of its bow shock. This 
situation cannot be "improved" if one takes into account that the 
true (ejection) velocity of HD 165319 (i.e. the peculiar velocity 
in the reference frame of the cluster) is higher than the "mea- 
sured" velocity (given in the last two columns of Table 1). Since 
the available proper motion measurements for this star are highly 
unreliable and the symmetry axis of its bow shock could deviate 
from the direction of stellar peculiar motion, one cannot how- 
ever exclude that either the true peculiar velocity of HD 165319 
is much higher than the "measured" one or that this star was 
ejected soon after the cluster formation. 

4. Discussion 

Our search for OB stars running away from the young open clus- 
ter NGC 661 1 has resulted in the discovery of three bow shocks 
generated by massive stars. The morphology of the bow shocks 
suggest that the associated stars were ejected from NGC 6611. 
All three stars are of O-type, which is consistent with the obser- 
vational fact that the percentage of O stars among runaways is 
higher than that of B stars (e.g. Stone TWTJ- 

It is natural to assume that the three stars were ejected when 
the central density in NGC 66 11 was sufficiently high for close 
dynamical encounters between the cluster members to be fre- 
quent. Whether or not the high star density in the cores of young 
clusters is primordial (i.e. intrinsic to the cluster formation pro- 
cess; e.g. Murray & Lin 1 19961 Clarke & Bonnell 2008) or is the 
result of the mass-segregation instability (e.g. Portegies Zwart 
et al. 119991 Giirkan et al. |2004| l remains unclear to date (e.g. 
Kroupa'2008). In the first case, the cluster could be collisional at 
birth, so that one can expect that the majority of runaway stars 
were produced at the very beginning of cluster evolution. In this 
case, the kinematic age of ejected stars should be comparable 
with the age of the cluster. In the second case, the beginning of 
the collisional stage is delayed for several Myr and the kinematic 
age of runaway stars is younger than the age of the parent clus- 
ter. Accurate proper motion measurements of runaway stars and 
their parent clusters are needed to distinguish between these two 
possibilities. 

The large separation between HD 165319 and NGC 6611 
and the low peculiar velocity of the star might imply that it 
was ejected soon after the cluster formation, supporting the first 
possibility. The proper motion measurements for this star are 
however unreliable (see Table 1), which leaves a possibility that 
HD 165319 was ejected recently, as suggested by the symmetry 
axis of its bow shock. The late production of runaway stars also 
follows from the young kinematic age of BD-I4° 5040 (the only 
one of the three stars with reliable proper motion measurements) 
and might be inferred from the morphology of the bow shock 
produced by star 1 . The available astrometric data, however, do 
not allow us to diff'erentiate between whether or not NGC 661 1 
experienced a burst of star ejection ~ 1 - 2 Myr ago or this 
process continued up until cluster formation. Indirect support of 
the latter possibility comes from the study of field O stars by 
Schilbach & Roser ( I2008I I. These authors retraced the orbit of 
the O-type star HD 157857 in the Galactic potential and found 



that this star was ejected from NGC 66 11 - 3.8 Myr ago with 
a peculiar velocity of ^ 114 km s"' . More robust proper motion 
measurements for star 1 and HD 165319 would help to solve the 
problem. 

By assuming that all three stars were ejected from NGC 661 1 
and taking into account that $ 20 per cent of runaway stars pro- 
duce bow shocks (see Sect. 1), one can infer that at least 15 
massive stars have been lost by the cluster during its lifetime. 
The true number of ejected OB stars could be far higher since 
high-velocity runaways would have already traveled beyond the 
region covered by the MSX survejQ, while the more numerous 
stars ejected at low velocities (^ 30kms"') have lower proba- 
bilities of their motion relative to the ambient medium being su- 
personic and therefore of being able to generate the bow shocks. 
We therefore believe that the percentage of stars producing bow 
shocks could be far lower than 20 per cent, based on the study 
of "classical" runaway stars. From this follows that the number 
of ejected OB stars could be comparable with or higher than the 
number still residing in NGC 661 1, which could explain the ap- 
parent deficit of OB stars in the cluster (see Sect. 2). 

In addition to the dynamical ejection, there are two other ef- 
fects that can reduce the size of the OB star population in the 
cluster: (i) supernova explosions of the most massive stars and 
(ii) merging of stars in the course of close binary-single and 
binary-binary star encounters. Supernova explosions in massive 
binaries could be accompanied by ejection of companion stars 
(e.g. Blaauw [1961| Leonard & Dewey IT993I I, which in addition 
reduce the population of massive stars in the cluster Supernova 
explosions begin to occur in clusters of age ^ 3 Myr. Flagey 
(i2007l l revealed a shell-like structure within the Eagle Nebula 
and suggested that this shell was a young (~ 3000 yr) super- 
nova remnant. If confirmed, this suggestion would imply that 
NGC 66 11 is at least as old as 3 Myr, which is consistent with 
the age of the cluster adopted in Sect. 2. For the age and mass 
of NGC 66 11 of ^ 4 Myr and ^ 2 x Mq, respectively, and 
assuming a Salpeter initial mass function, one finds that ^ 10 
stars have exploded as supernovae since the cluster formation. 
Correspondingly, the number of massive runaways produced 
by binary disruptions following supernova explosions cannot 
exceed this number. Merging of stars is a by-product of star 
ejection by means of binary-single and binary-binary encoun- 
ters (e.g. Leonard 1995). Four early-type O stars in NGC 6611 
and the star producing the bow shock 2 (BD- 14° 5040) appear 
younger than the main population of OB stars in the cluster and 
therefore could be the rejuvenated merged stars (e.g. Leonard 
1995; Portegies Zwart et al. '19991 Gvaramadze & Bomans 
2008). Therefore, we suggest that at least 5 OB stars were "lost" 
by the cluster via merging. 

We note that HD 165319 was considered by de Wit et al. 
(120051 1 to be one of "the best examples for isolated Galactic high- 
mass star formation", de Wit et al. ( 120041 1200 5) searched for 
parent star clusters for 43 O-type field stars to confirm whether 
or not these massive stars originated in a clustered mode of star 
formation. They found that 1 1 stars from their sample cannot 
be associated with a nearby cluster or star-forming region and 
suggested that ~ 4 + 2 per cent of all field O-type stars could 
be formed in isolation (cf. Chu & Gruendl 2008). If confirmed, 
this result could have important implications for understanding 
massive star formation. It would either imply the existence of a 
special mode of isolated star formation (e.g. Yorke & Sonnhalter 
120021 Li et al. |2003i l or allow the possibility that low mass (< 



' A possible example of such a high-velocity runaway is the afore- 
mentioned star HD 157857, located at ~ 13° from NGC 661 1. 
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Fig. 5. A 0°5 X 0°.5 SHS image of RCW 158, again oriented in 
Galactic coordinates (increasing 11 to left, and increasing b up- 
ward). 



100 Mq) clusters are able to form O stars (e.g. Elmegreen |2006l 
Parker & Goodwin 2007 ). 

de Wit et al. (2005) restricted the search for parent clus- 
ters within a circle of radius 65 pc; this figure comes from 
the average peculiar radial velocity of field O stars, which is 
6.4kms~' (Gies 11987!), multiplied by their average lifetime of 
10 Myr (1 kms"' x 1 Myr ^ 1 pc). It is possible therefore 
that at least some of the 1 1 stars, suggested by de Wit et al. 
(120051 ) to be formed in isolation, were instead ejected from clus- 
ters more distant than 65 pc. The discovery of the bow shock 
around HD 165319 and identification of its possible parent clus- 
ter (separated by ^ 105 pc) lends support to this suggestion. We 
note that some of the remaining 10 stars could be rejuvenated 
stars (blue stragglers) formed via merging of less massive stars 
in binaries ejected from clusters located far beyond the circle of 
radius of 65 pc. Study of field O-type stars by Schilbach & Roser 
( 1200 8 ) suggests that 6 of these 10 stars were ejected from open 
star clusters and were therefore formed in the clustered way. In 
some cases, however, the times-of-flight derived by Schilbach 
& Roser (120081 ) exceed the lifetimes of the stars. The proposed 
identification of the parent clusters would therefore be incorrect 
if these stars are not blue stragglers. 

We searched for bow shocks around the remaining 10 O- 
type stars formed in 'isolation' using the MSX survey but none 
were found (note that the survey covers only 8 of 10 stars). This 
non-detection could be due to the limited and inhomogeneous 
sensitivity of the survey (see Price et al. 120011) . It is also possible 
that some of these stars have not produced bow shocks simply 
because their peculiar velocities are less than the sound speed 
in the (hot) ambient interstellar medium (see Sect. 1). We note 
that the peculiar velocities relative to the ambient medium (and 
therefore ability to generate bow shocks) of some runaway stars 
could be affected by the cluster peculiar motion. A (massive) star 
ejected in the direction opposite to the cluster's motion will have 
a reduced peculiar velocity relative to the local standard of rest, 
while its true (ejection) velocity with respect to the cluster could 
be sufficiently high to classify the star as a "classical" runaway. 
This star has not only a lower probability of producing a bow 
shock, but is also more distant from the parent cluster than one 
can infer from its apparent proper motion. Since most ejected 
stars have velocities of several tens km s"' (i.e. comparable with 
the escape velocity from the parent cluster) and the peculiar ve- 
locities of clusters could be as high as ~ 10 - 20kms"' (de 
Zeeuw et al. 1999; see also Sect. 2), the above effect could be 
significant and might cause an excess of bow-shock-producing 
stars ahead of the moving cluster 

In Sect. 3. 1, we mentioned that the bow shock 3 is part of the 
H II region RCW158. de Wit et al. ( |2005| ) consider this nebula as 
the tracer of isolated high mass star formation. The almost cen- 
tral location of HD 165319 in RCW 158 (Fig. 5) suggests, how- 
ever, that RCW 158 might be a Stromgren zone created by the 
ultraviolet emission of the star. We now check this possibility. 
Assuming that the gas in RCW 158 is fully ionized and its tem- 
perature is 10"* K, one can obtain an estimate of the gas number 
density (e.g. Lequeux l2005l ) 



n ^ 22 cm 



-3/2 



S(0) 



6.88 pc / \4.90 X 10*^ photons s" 



1/2 



where Rs is the Stromgren radius (derived from the angular ra- 
dius of RCW 158 of ~ 12' and t/ = 2 kpc) and 5(0) is the to- 
tal ionizing-photon luminosity of an 09.5Iab star (taken from 
Martins et al. 2005 ). This value should be equal to that derived 
from the ram pressure balance between the stellar wind and the 
ambient medium: 



n ^ 17 cm 



M 



5.3x lO-''M0yr- 



l)(21OOkms-l) 



""(osi^) (50kms-i) ' 

where M and Voo are the wind mass-loss rate and terminal ve- 
locity of an 09.5Iab star (taken from Mokiem et al. 120071 ), Rq 
is the stand-off radiufl and v, is the velocity of the star rela- 
tive to the ambient medium, assumed to be equal to 50kms"' 
(see Sect. 3.3). The good agreement between the two estimates 
supports our suggestion that RCW 158 is the Stromgren zone 
of HD 165319. The stellar parameters used in Eq. 2 could how- 
ever still vary within the classification [see e.g. the values given 
by Markova et al. ( 120041 ) for the same spectral and luminosity 
class]. Using different possible values for the stellar parameters 
and allowing the velocity to be between 30 km s"' and 60 km s"' 
in Eq. 2, we estimate the approximate lower and upper bounds 
of the density to be 10 cm"^ and 30 cm"-', respectively. This 
density range is again in good agreement with the density esti- 
mate from the Stromgren formula (Eq. 1). Further support comes 
from the observation that RCW 158 is partially surrounded by 
8.3 fim filamentary structures, which are typical of photodisso- 
ciation regions (Hollenbach & Tielens'1997). We therefore con- 
clude that the ionizing photons emitted by HD 165319 have cre- 
ated a Stromgren sphere in the dense ambient medium (visible 
as RCW 158), while the interaction between the wind of this su- 
personically moving star and the ambient gas generates a bow 
shock ahead of the star and evacuates a cavity behind (see Fig. 5; 
cf. Raga et al. .l997j . 

5. Conclusions 

With the detection of three bow-shock-producing stars in the 
vicinity of NGC 66 1 1 , we have found strong evidence that the 
proposed loss of massive stars, due to dynamical processes in 
the early evolution of young clusters, is indeed operating. While 
existing astrometric data for the three stars do not allow us to 
determine precisely the timing of the ejection or study the inter- 
action of the runaway stars with the ISM in detail, the morphol- 
ogy of the bow shocks and the proper motion of the cluster itself 
allows us to argue that these stars originated in NGC 66 1 1 . 

Clearly, our finding opens up a number of possibilities for 
follow-up investigations of the NGC 66 11 bow shocks and the 
ejected stars. More accurate proper motion measurements will 
allow us to check if the stars were ejected continuously or dur- 
ing a short period of cluster evolution. The most important point 
nevertheless is that the idea proposed by Pflamm-Altenburg & 
Kroupa (2006) is supported by direct observation, at least in the 
case of NGC 6611. The results for several other young clusters 
will be presented in forthcoming papers. 
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(1) 



- We neglected here a geometric factor of order unity, taking into ac- 
count the inclination of the bow shock due to the peculiar radial velocity 
of the star of =^ 20.5 km s"' [derived from the measured radial velocity 
of ^ 25.4 km S-' (Kharchenko et al. l2007] )1. 
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